David  Taylor  Research  Center 

Bethesda,  Maryland  20084-5000 


Carderock  Div/SHD- 1387-02  December  1992 

Ship  Hydromechanics  Department 


Z 

g 

< 


O 

CJ 

Q 

W 

W 

Ph 

CO 


a 

a 

M 

a 

a 

o 

r 


2 

O 

M 

H 

U 


Q 

W 

a 

p* 


c 

M 

E- 

U 

a 

c 

w 

Q 

Eh 

CO 


a 

a 

a 

EH 


CM 

O 


r-~ 

00 

co 


Q 

a 

CO 


> 

•H 

o 

a 

a 

O 

<u 

■a 

Vi 

a 


THRUST  DEDUCTION  PREDICTION  FOR  HIGH  SPEED 
COMBATANT  SHIP 


By 

Yoon-HoKim 


DTIC 


ELECTE 
JUNl  5  1993 


Approved  for  Public  Release,  Distribution  Unlimited 


93-13296  0 


Bs  6  14 


REPORT  DOCUMENTATION  PAGE 

la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

lb  RESTRICTIVE  MARKINGS 

2a  SECURITY  CLASSIFICATION  AUTHORITY 

3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

2b  DECLASSIFICATION /OOWNGRAr'NG  SCHEDULE 

Approved  for  Public  Release; 

Distribution  Unlimited 

4  PERFORMING  OPOANI2ATION  REPORT  NUMBER(S) 

Carderock  Div/SHD  1387-02 

S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6a  NAME  OF  PERFORMING  ORGAN I2ATION 
David  Taylor  Model  Basin 
Carderock  Division  NSWC 

6b  OFFICE  SYMBOL 
(If  applicabia) 

7a.  NAME  OF  MONITORING  ORGANIZATION 

6c  ADDRESS  (Gty.  Start,  and  IIP  Coda) 

Bethesda,  MD  20084-5000 

7b.  ADDRESS  (City.  Start,  and  HP  Coda) 

6a  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 

Naval  Sea  Systems  Command 

8b.  OFFICE  SYMBOL 
(If  appUcabta)- 

052 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

Washington,  D.C.  20362-5101 


tO  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT 

ELEMENT  NO.  NO 

63573N 


tl  HUE  (includt  Security  Clauitication) 

Thrust  Deduction  Prediction  for  High  Speed  Combatant  Ship 


1 3£inal  0r  K*°*T 


13b  TIME  COVERED 
FROM  _______  TO 


um^c^w’Mont,KDs,> 


COSATI  COOES 


FIELD  I  GROUP  SUB-GROUP 


II.  SUBJECT  TERMS  IContmt*  on  rtvana  if  nma tty  and  Kttrnify  fey  fetorf  number) 
Resistance  Actuator  Disk 

Powering  Thrust 

Panel  Method _ _ . 


19  ABSTRACT  (Continue  on  reverie  if  nttauaty  and  tdtnttfy  fey  feJodt  number ) 

The  free  surface  effects  on  the  huli-propeiler  interaction  characteristic  of  thrust  deduction  have  been  studied  for  a 
high  speed  combatant  ship  hull  form.  Through  mathematical  modelling  and  a  series  of  numerical  computations, 
we  attempt  to  elucidate  the  mechanism  of  this  component  of  hull-propeller  interaction  from  the  free  surface. 
SWIFT,  a  linearized  free  surface  potential  flow  solver  using  a  higher-order  panel  method,  has  been  extended  to 
study  the  interaction  between  the  double  hull  and  propeller  as  well  as  the  hull-propeller  interaction  under  free 
surface  waves.  The  propeller  is  simulated  using  an  actuator  disk.  Thrust  loading  coefficients  from  propeller  open- 
water  tests  provide  die  sink  strength  for  the  disk.  This  ship  has  inclined  propeller  shafts  and  the  measured  wakes 
at  die  propeller  plane  clearly  indicate  that  the  propellers  are  operating  outside  the  turbulent  boundary  layer  at  the 
ship  stem.  The  viscous  effects  arc  not  treated  here.  The  pressure  on  the  hull  is  integrated  to  obtain  the  resistance 
for  the  double  body  flow  with/without  propellers,  and  also  for  the  body  and  free  surface  flow  with/without 
propellers.  The  computed  thrust  deduction  fractions  show  good  agreement  with  the  experimental  data.  For  this 
particular  ship,  the  free  surface  wave  effect  on  thrust  deduction  is  not  small  and  lies  between  25%  to  30%  of  the 
total,  depending  upon  speeds. 


20  DISTRIBUTION/ AVAILABILITY  OF  ABSTRACT 
G UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT  □  OTIC  USERS 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Yoon-HO  Kim 


DD  FORM  f 473, 84  MAR  43  apa  turnon  may  bt  uwd  until  amautttU 

All  othtr  td'tiont  are  obioltta 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b  TELEPHONE  (kKlodt  PraaCodt)  22c  OFFICE  SYMBOL 

(301)227-1594  code  1522 _ 


ii  amautttU  security  classification  of  this  page 


CONTENTS 


NOMENCLATURE 

ABSTRACT 

ADMINISTRATIVE  INFORMATION 
INTRODUCTION 
NUMERICAL  APPROACH 


Ship,  Hull 

Propeller 

Free-Surface  Wave 

COMPONENTS  OF  THRUST  DEDUCTION 


Double  Model  Potential  Fraction 
Free  Surface  Wave  Fraction 


DESCRIPTION  OF  MODEL 
RESULTS  AND  DISCUSSION 
CONCLUSIONS  AND  RECOMMENDATIONS 
ACKNOWLEDGMENTS 
REFERENCES 


Accesion  For 


NTIS  CRA&I 
DTIC  TAB 
Unannounced 
Justification 


By  . . 

Distribution/ 


J 


LJ 


Availability  Codes 

Dist 

&L 

Avail  a 
Spe 

nd  /  or 
cial 

iii 


FIGURES 


1 .  Coordinate  System  1 1 

2.  Typical  Panel  Arrangement  for  a  Ship  and  Free  Surface  12 

3.  Propeller  Location  and  Dimension  1 3 

4.  Ship  Hull  Panel  Network  14 

5.  Comparisons  of  Cross  Flow  Velocity  Vectors  with  and  without  an  Actuator  Disk  15 

6.  Comparison  of  Wave-Making  Resistance  Coefficients  Predicted  by  SWIFT  and 

from  Model  Experiments  17 

TABLES 

1.  Ship  Particulars  18 

2.  Mean  Axial  and  Total  Velocity  Increment  due  to  Actuator  Disk  at  Fn=0.41 4  1 9 

3.  Potential  Fraction  of  Thrust  Deduction  20 

4.  Free  Surface  Wave  Fraction  of  Thrust  Deduction  2 1 

5.  Thrust  Deduction  Comparison  between  Experiment  and  Computation  22 


iv 


NOMENCLATURE 


A<> 

CD 

CP 

CR 

^TH 

D 

g 

h 

q 

R 

Rp 

Rh 

s 

T 

t 

tp,tv>tw 


U 

v 

V 

VA 


we 

Wp,Wv,Ww 

wj 

x,y,z 


Propeller  disk  area,  Aq  =  tc(Rp2  -  Rh2) 

Drag  coefficient,  CD  =  D/(^pU2A0) 

Pressure  coefficient,  C?  =  Pressure/(-jpU2) 

Resistance  coefficient,  CR  =  Resistance/ (^pU2S) 

Thrust  loadiig  coefficient,  =  T/(|pA0U2) 

Drag  augmentation,  D  =  T  -  R 
Gravitational  acceleration 

Transom  stem  depth  measured  from  undisturbed  free  surface  to  bottom  edge 

Sink  strength  at  an  actuator  disk 

Barehull  Resistance 

Propeller  tip  radius 

Propeller  hub  radius 

Ship  wetted  surface 

Propeller  thrust 

Thrust  deduction  fraction 

Thrust  deduction  fraction  due  to  double  body  potential,  viscous,  and  free- 

suxface  wave,  respectively 

Ship  speed 

Pertur  bed  velocity 

Total  velocity,  V=U+v 

Speed  of  advance  of  the  propeller 

Wake 

Effective  wake 

Wake  due  to  double  body  potential,  viscous,  free  surface  wave,  respectively 
Taylor  wake  fraction,  determined  by  measurement  in  a  powering  experiment 
Cartesian  coordinate  system  fixed  on  the  ship 


d>  Total  velocity  potential 

0  Perturbed  velocity  potential 

^db^h^w  Velocity  potential  for  double  body,  ship,  and  free  surface,  respectively 
<Hp  Velocity  potential  due  to  propeller 

p  Fluid  density 

o(^,T])  Source  strength  distribution  across  a  panel 

X  Local  coordinate  system 


ABSTRACT 


The  free  surface  effects  on  the  hull-propeller  interaction  characteristic  of  thrust  deduction 
have  been  studied  for  a  high  speed  combatant  ship  hull  form.  Through  mathematical 
modelling  and  a  series  cf  numerical  computations,  we  attempt  to  elucidate  the  mechanism  of 
this  component  of  hull-propeller  interaction  from  the  free  surface.  SWIFT,  a  linearized  free 
surface  potential  flow  solver  using  a  higher-order  panel  method,  has  been  extended  to  study  the 
interaction  between  the  double  hull  and  propeller  as  well  as  the  hull-propeller  interaction  under 
free  surface  waves.  The  propeller  is  simulated  using  an  actuator  disk.  Thrust  loading  coeffi¬ 
cients  from  propeller  open-water  tests  provide  the  sink  strength  for  the  disk.  This  ship  has 
inclined  propeller  shafts  and  the  measured  wakes  at  the  propeller  plane  clearly  indicate  that  the 
propellers  are  operating  outside  the  turbulent  boundary  layer  at  the  ship  stem.  The  viscous 
effects  are  not  treated  here.  Tne  pressure  on  the  hull  is  integrated  to  obtain  the  resistance  for 
the  double  body  flow  with/without  propellers,  and  also  for  the  body  and  free  surface  flow 
with/without  propellers.  The  computed  thrust  deduction  fractions  show  good  agreement  with 
the  experimental  data.  For  this  particular  ship,  the  free  surface  wave  effect  on  thrust  deduction 
is  not  small  and  lies  between  25%  to  30%  of  the  total,  depending  upon  speeds. 
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INTRODUCTION 

When  two  bodies  in  close  proximity  to  each  other  move  in  a  fluid  there  are  interaction 
effects  which  cause  additional  forces  and  moments  to  be  experienced  by  each  body.  The  additional 
force  arising  from  propeller-hull  interaction  is  well  known  to  naval  architects.  A  propeller  operating 
at  the  ship  stem  accelerates  the  local  flow  over  the  stem  of  the  ship  and  results  in  changing 
pressure  distribution  which  has  effect  on  the  hull  drag.  This  change  in  drag  of  the  hull  is  called 
drag  augmentation  and  is  expressed  in  terms  of  the  thrust  deduction.  The  flow  acceleration  also 
increases  the  wall  shear  stress,  and  hence,  the  frictional  drag.  The  flow  interactions  have  important 
consequences  to  the  engineering  estimates  of  speed  and  power  for  the  ships.  For  instance,  flow- 
related  features  such  as  propeller  proximity  to  the  hull,  local  hull  shape,  overall  body  flow,  and  the 
free  surface  influence  are  all  involved  with  determining  net  thrust  required  and  understanding  the 
propulsor-induced  drag  penalty.  Details  and  order  of  magnitude  estimates  are  needed  for 
preliminary  design  studies.  Our  efforts  are  concentrated  on  clarifying  the  quantitative  role  of 
wavemakir.g  due  to  the  presence  of  the  free  surface  in  the  phenomenon  of  hull-propulsor 
interaction  and  consequently  its  contribution  to  the  hydrodynamic  propulsive  efficiency  of  the  hull- 
propulsor  system. 

The  interaction  may  be  conveniently  studied  in  terms  of  propulsion  factors,  including  wake 
and  thrust  deduction  fraction.  The  wake  w  caused  by  the  presence  of  the  hull  and  the  free  surface 
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is  a  simple  measure  of  the  change  in  propulsor  inflow  as  compared  to  an  equivalent  open -water 
condition,  and  is  defined  as  follows: 


(1) 


where  U  is  the  ship  speed  and  Va  is  the  speed  of  advance  of  the  propeller. 

The  thrust  deduction  fraction  t  is  an  indirect  expression  of  the  fact  that  the  force  of 
resistance  acting  on  the  hull  is  modified  as  a  result  of  propeller  action  and  is  defined  as  follows: 


T-R 

t  =  - 

T 


(2) 


where  R  is  the  barehuii  resistance  and  T  is  the  propeller  thrust.  Following  Dickmann's  pioneering 
work[l],  the  hull-propulsor  interaction  has  been  customarily  considered  as  a  superposition  of  three 
basic  effects:  zero-Froude-number  potential  effects,  viscous  effects,  and  wave  effects  due  to  the 
presence  of  the  free  surface.  Using  standard  notations,  wake  and  thrust  deduction  fractions  are 
expressed  as 


w=  Wp+  Wv  +  ww 

(3) 

t  =  tp  +  tv  +  tw 

(4) 

where  the  subscripts  p,  v,  and  w  denote  potential,  viscous,  and  free  surface  wave,  respectively. 

Since  Dickmannfl],  considerable  effort  has  been  made  to  study  the  potential  and  viscous 
effects  in  hull  propeller  interaction  (see  the  comprehensive  bibliography  by  Nowacki  and 
Sharma[2]).  Beveridge^]  was  the  first  to  apply  a  panel  method  to  represent  the  practical  three- 
dimensional  body  and  to  distribute  Rankine  type  sources  across  each  panel  to  calculate  the  flow. 
He  used  a  constant  strength  sink  disk  for  the  propeller  and  later  introduced  a  sink  disk  with  radially 
varying  strength  derived  from  propeller  lifting-line  theory  for  contrarotating  propellers.  He 
computed  thrust  deduction  using  Lagally's  theorem  and  the  results  showed  good  agreement  with 
the  experimental  data.  Stem  et  al[4]  described  a  comprehensive  viscous-flow  method  for  the 
computation  of  propeln  i-hull  interaction  in  which  a  numerical  method  for  calculating  the  viscous 
flow  over  the  stem  and  in  the  wake  of  a  ship  is  coupled  with  a  propeller-performance  in  an 
interactive  and  iterative  manner  to  predict  the  combined  flow  field.  In  their  computation  the  free 
surface  effect  was  not  included  and  furthermore  a  propeller  is  represented  approximately  by  its 
induced  velocity  effects.  Therefore,  for  specified  propeller  geometry,  operating  condition,  and 
nonuniform  inflow-velocity  distribution,  the  propeller  induced  effects  are  assumed  known. 

Nowacki-Sharma[2]  investigated  the  wave  effects  due  to  the  presence  of  the  free  surface 
both  experimentally  and  analytically.  Thin  ship  theory  was  used  to  compute  wave-making 
resistance.  The  propeller  induced  flow  was  approximated  by  distributing  Havelock  source  on  the 
actuator  disk.  They  showed  that  the  free  surface  wave  component  can  be  dominant  in  the  wake  and 
quite  significant  in  the  thrust  deduction  at  some  Froude  numbers.  Yamazaki-Nakatake[5] 
conducted  experiments  using  the  similar  mathematical  hull  of  Nowacki-Sharma.  They  used  a 
different  model  propeller,  but  kept  the  same  propeller-diameter  to  the  ship-draft  ratio.  With  careful 
and  systematic  series  of  experiments  they  observed  that  the  thrust  deduction  fraction  generally 
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varies  with  the  extent  of  wave  making.  All  the  studies  showed  that  the  wave  making  effect  on  the 
thrust  deduction  is  not  small,  even  at  high  speed. 

The  wake  fraction  for  high  speed  combatant  ships  is  known  experimentally  to  be  small, 
usually  between  ±0.02,  while  the  wake  fraction  for  single  screw  commercial  ships  varies  from  0.2 
for  low  block  to  0.5  for  high  block  coefficient  ship.  The  ship  hull  we  consider  in  this  study  is  a 
high  speed,  transom  stem  ship.  The  measured  wake  fraction  in  the  propeller  plane  for  this  model 
ranges  from  0.01  to  0.04  depending  upon  the  speeds.  A  panel  method  is  employed  for  this  study. 
The  zero-Froude-number  potential  effects  are  computed  first.  Next  when  we  consider  the  free 
surface  wave  effects,  the  body  boundary  condition  has  been  modified  at  the  transom  stem.  In  the 
following,  we  briefly  outline  the  numerical  approach  and  the  computed  results  are  discussed  along 
with  the  measurements. 


NUMERICAL  APPROACH 

Suppose  that  a  ship  moves  in  the  positive  x-dinection  with  constant  forward  speed  U  in 
calm  water.  A  Cartesian  coordinate  system  xj  =  (x,y,z)  translates  with  the  ship  as  shown  in  Fig.  1. 
The  z=0  plane  is  taken  as  the  undisturbed  free  surface,  the  positive  x-axis  in  the  direction  of  the 
ship's  forward  velocity,  and  the  positive  z-axis  upward.  The  boundary-value  problem  will  be 
expressed  relative  to  this  moving  coordinate  system  with  the  flow  at  infinity  consisting  of  an 
uniform  stream.  The  fluid  is  assumed  to  be  inviscid  and  incompressible  and  its  motion  is 
iirotational  such  that  the  velocity  field  of  the  fluid  can  be  defined  as 

v(x,y,z)  =  V<j>(x,y,z)  (5) 

where  <|>(x,y,z)  is  the  velocity  potential  and  satisfies  the  Laplace  equation 

=  0  <6) 


in  the  fluid  domain  D  and  the  body  boundary  condition 

V<j)-n  =  U-n  (7) 

on  the  body  surface  Sb  where  n  =  (nx,  ny,  nz)  denotes  the  outward  unit  normal  vector  on  the 
boundary.  Furthermore  the  disturbance  due  to  free  surface  is  assumed  to  be  sufficiently  small  that 
the  nonlinear  free  surface  condition  may  be  linearized. 

+  g¥,  =  0 .  at  z=0  (8) 

Energy  conservation  requires  that  velocity  potential  approach  the  uniform  onset  flow  potential  and 
that  there  be  no  waves  far  upstream  of  the  ship,  and  that  waves  always  travel  downstream. 

There  are  two  common  approaches  to  solve  this  boundary  value  problem :  Havelock  source 
and  Rankine  type  singularity  distribution  methods.  Havelock  source  has  been  developed  to  satisfy 
the  linearized  free  surface  and  radiation  conditions.  But,  because  of  its  complex  kernel  function,  its 
application  to  practical  hull  configurations  is  much  more  difficult  and  cumbersome.  Rankine  type 
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singularity,  on  the  other  hand,  has  simple  kernel  function  and  is  easier  for  handling  the  details  of 
complex  hull  geometry.  The  geometry  of  ship  stem  fitted  with  propulsion  system  is  particularly 
complicated  and  a  boundary  element  method  may  be  the  only  method  capable  to  handle  the  details 
of  its  complex  configuration  with  sufficient  accuracy.  A  well-known  drawback  of  this  approach  is 
the  necessity  of  proper  numerical  implementation  of  the  radiation  condition  for  free  surface  waves. 
Dawson[6]  successfully  introduced  an  upstream  finite  difference  operator  to  Ranldne  method  when 
implementing  the  free  surface  condition  in  order  to  satisfy  the  radiation  condition. 

The  method  adopted  here  is  a  boundary  element  method.  The  ship  hull,  propeller,  and  free 
surface  are  discretized  into  many  small  quadrilaterals.  Across  each  qaadrilateral  the  Rankine 
source(lfr)  is  distributed.  The  source  strength  o  at  a  point  (£,T|,£)  is  assumed  to  vaiy  linearly 
across  the  quadrilateral 

a($,Ti)  =  a0+a^  +  a„ti.  (9) 

The  coefficients  in  Eq.  (9)  are  determined  by  using  a  method  of  weighted  squares  over  the  panel 
and  up  to  eight  of  its  neighboring  panels.  This  method  requires  that  the  form  of  Eq.  (9)  give 
exactly  the  singularity  value  at  its  centroid  and  approximate  values  in  the  least  square  sense  at 
centroid  of  the  neighboring  panels.  We  used  curved  panels  to  patch  the  ship  ana  free  surface.  Each 
quadrilateral  may  be  represented  in  the  approximate  form 

Cfeti)  -  (10) 

where  (^,T],Q  are  orthogonal  coordinates  local  to  S.  The  six  coefficients  in  Eq.(10)  are  determined 
by  requiring  that  the  approximate  surface  given  by  Eq.  (10)  pass  through  its  four  comer  points 
exactly  and  through  its  neighboring  points  approximately  in  a  least  square  sense.  The  boundary- 
value  problem  formulated  above  then  reduces  to  a  determination  of  an  unknown  singularity 
distribution  over  the  boundary  surface  of  the  fluid  domain. 

We  linearized  the  boundary  value  problem  and  hence  the  velocity  potential  (j>  can  be 
assumed  linear  superposition  of  each  disturbance.  In  case  of  a  double  body  model  without 
propeller  0  =  0db,  a  double  body  with  propeller  0  =  0db+0p,  a  towed  ship  without  propeller  <)>  = 
<t>h+<J>w»  and  a  self-propelling  ship  0  =  0h+0w+0p-  In  the  following,  each  velocity  potential  will 
be  discussed  first  and  components  of  thrust  deduction  fraction  tp  and  tw  will  be  considered  next 


Ship  Hull 

The  velocity  potential  induced  by  the  Rankine  source  distribution  on  the  hull  surface  may 
be  given  by 


4>h 


—  f  2ds 
4  k  S.  r 


or  for  a  double  body  model  case 


(11) 

(11)' 
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where  r  =  ( (x-£)2+(y-Tl)2+(z-02  )l/2,  j>  =  ( (x-£)2+(y-Tl)2+(z+C)2  )l/2,  c  the  unknown  source 
strength,  and  Sb  the  wetted  hull  surface.  If  we  discretize  the  hull  into  Nb  number  of  panels,  then 
the  velocity  potential  can  be  expressed  as  follows: 


“  ■  -il  l  7* 


•■--iKHW) 


The  unknown  singularity  strength  distribution  c  is  determined  by  satisfying  the  body  boundary 
condition  (7). 

Propeller 

Following  Dickmannfl],  the  propeller  is  modelled  by  a  sin-flow  actuator  disk.  Similar  to 
Eq.(l  1)  the  propeller  induced  velocity  potential  may  be  expressed  as 


= 


-Lj s5ds  „  iffj  1^ 

JspT  4it  1rfvJs<r  J 


where  Np  is  the  total  number  of  panels  used  to  represent  the  propeller.  In  case  of  a  double  body 
model,  we  use  z=0  as  a  plane  of  symmetry  and  have  a  similar  expression  of  Eqs.  (1 1)’  and  (12)’. 
Unknown  sink  strength  q  is  assumed  constant.  Propeller  loading  is  incorporated  when  determining 
unknown  sink  strength  distribution  q  across  the  disk  as  follows: 

q  =  (l-wE)  -  (14) 

where  we  is  the  measured  effective  wake  fraction.  The  thrust  loading  coefficient  Cth  is 
customarily  given  as 

^  =  IpM?  (15) 

where  p  is  the  water  density,  A0  =  7t(  Rp2-Rh2)  the  propeller  disk  area  in  which  Rp  and  Rh  are  the 
propeller  tip  and  the  propeller  hub  radius,  respectively.  The  drag  augmentation  D  is  the  difference 
between  the  thrust  and  the  resistance,  i.e.,  D  =  T-R  and  its  coefficient  is  defined  by 


W,  -  1 


ipA0U 
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From  Eqs.  (2),  (15),  and  (16)  the  thrust  deduction  fraction  t  can  be  expressed  as  follows: 


T-R  D  _  CD 
T  =  T  C™ 


Free  Surface  Wave 

The  velocity  patential  due  to  the  presence  of  the  free  surface  has  the  similar  expression  of 
Eqs.(ll)  and  (12),  i.e., 


(18) 


Here  infinite  free  surface  domain  is  truncated  into  the  finite  domain  Sw  for  computational  reason, 
and  Nw  is  the  total  number  of  panels  used  to  represent  the  truncated  free  surface  domain. 

DDG  51  hull  has  transom  stem.  During  the  model  experiment  we  observed  that  at  ship 
speed  greater  than  20  knots  the  flow  clears  the  stem  completely  and  the  cross  section  becomes  dry. 
In  order  to  model  this  physical  phenomena  numerically,  we  impose  two  conditions  at  the  stem. 
First,  to  satisfy  the  dry  condition  at  the  cross  section  of  the  transom  stem,  we  do  not  put  any 
panels  to  represent  the  cross  section  and  leave  it  empty,  in  other  words,  the  ship  hull  is  not  closed 
at  the  stem.  Secondly,  at  the  bottom  edge  of  the  transom  stem,  the  body  boundary  condition  of 
Eq.(7)  is  replaced  by  the  requirement  that  the  flow  pass  by  each  bottom  edge  panel  tangentially. 
The  magnitude  of  the  tangential  velocity  is  assumed  equal  to  square  root  of  the  hydrostatic  term 
2gh  of  Bernoulli  equation 


P  2 


+ 


gh 


Po 

P 


(19) 


in  which  p0  is  the  atmospheric  pressure,  V  =  U  +  v  t  g  js  the  gravitational  acceleration,  and  h  is  the 
depth  of  the  transom  stem  measured  from  the  undisturbed  free  surface  to  the  bottom  edge.  The 
angle  of  tangential  velocity  is  obtained  by  extrapolating  values  from  neighboring  panels  at  the 
bottom  edge  of  transom  stem. 

A  computer  code  named  SWIFT(£urface  Wave  Inviscid  Flow  Theory),  using  a  higher 
order  panel  method,  has  been  developed  by  Kim  et  al.[7]  to  solve  the  surface  ship  problem  with 
the  linearized  free  surface  condition.  Numerous  computed  results  demonstrate  its  capability  to 
accurately  predict  the  flow  field,  wave  profiles,  and  wave-making  resistance  for  cruiser  stem  af 
well  as  transom  stem  type  ships.  The  SWIFT  computer  code  has  been  used  extensively  for  this 
study. 
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COMPONENTS  OF  THRUST  DEDUCTION 


Double  Model  Potential  Fraction 

The  velocity  potential  for  the  double  body  model  $  =  <(>db  of  (1 1)'  satisfies  the  Laplace 
equation  (6)  automatically  and  unknown  singularity  strength  o  at  each  panel  can  be  determined  by 
satisfying  the  body  boundary  condition  (7).  Once  the  singularity  distributions  are  determined  as  the 
solution  of  the  boundary- value  problem,  the  velocity  potential  can  be  obtained  using  Eq.(l  1)'  at 
any  point  in  the  fluid  domain  and  the  velocity  field  is  determined  from  Eq.(5).  The  hydrodynamic 
pressure  coefficient  is  given  by  Bernoulli  equation(19): 

P  V2 

CP  =  — ^  =  1~ 

ipu2  U2  (20) 

and  the  resistance  coefficient  Cr  can  be  computed  by  integrating  the  pressure  over  the  wetted  hull 
surface 


(21) 


In  case  of  a  double  body  model  fitted  with  propeller(s),  the  total  velocity  potential  <i> 
consists  of  two  terms,  <t>db  and  <|>p.  The  presence  of  propeller(s)  must  be  taken  into  account  when 
satisfying  the  body  condition  (7)  to  determine  the  strength  of  singularities  distributed  on  the  hull. 
The  body  boundary  condition  accordingly  becomes 

B  =  -Un,  -  V*p  n  on  Sb  (22) 


Again  after  determining  the  singularity  strength  distribution  a  on  the  hull,  we  compute  the  velocity 
distribution  on  the  body  from  Eq.(S),  the  pressure  distribution  from  Eq.(20),  and  finally  the 
resistance  coefficient  from  Eq.(21).  The  potential  fraction  tp  is  then  obtained  from 


t 


p 


(CrX/Ptop 


w/oProp 


(“TH 


(23) 


where  (C?)w/ft0p  and  (C*)w/oProp  are  the  double  body  resistance  coefficients  with  and  without 
operating  propellers),  respectively. 

Free  Surface  Wave  Fraction 


The  velocity  potential  for  a  towed  ship  may  be  expressed  as  (j>  =  <j»h  +  <J>W-  Unknown  o's 
are  determined  by  satisfying  both  body  boundary  condition  (7)  and  free  surface  condition  (8). 
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Figure  2  shows  the  typical  panel  network  used  for  the  computation.  For  a  self-propelling  ship,  the 
total  velocity  potential  <)>  consists  of  three  components:  4»j,»  4*w>  and  4>p.  Here  <J>P  is  assumed 
known.  To  determine  the  singularity  strength  distribution  o  the  boundary  conditions  have  to 
change  accordingly,  and  the  body  boundary  condition  becomes 

V(<J>h  +  <t>„)-n  =  -Unx  -  V(J>  -  n  onSe 

and  the  free  surface  condition  becomes 

U2(4>h+<U,*  +  g(4>h+<U  =  -U2(4>p)„  -  gC^p),  on  z=0 


After  determining  a,  computations  of  velocity  potential  and  velocity  field  as  well  as  pressure  and 
resistance  are  straightforward,  and  the  desired  wave  fraction  tw  is  obtained  as  follows: 


t_  = 


(Cr  )»/Prop 


(Ct') 


R  'wloProp 


"TH 


(24) 


where  (C^*)wlPtop  and  (C^m)vloftop  are  the  resistance  coefficients  when  a  ship  is  moving  in  calm 
water  with  and  without  operating  propellers),  respectively. 


DESCRIPTION  OF  MODEL 


Table  1  shows  the  particulars  of  the  ship  hull.  All  the  quantities  are  nondimensionalized.  A 
dome  is  installed  at  the  bow  and  the  stem  type  is  transom.  The  propeller  shaft  is  inclined  4.5 
degree  with  respect  to  the  keel  line.  Figure  3  is  a  schematic  drawing  of  the  propeller  location,  tip 
radius,  and  the  hub  radius.  All  the  model  experiments  were  conducted  in  the  deep  water  towing 
tank  of  David  Taylor  Model  Basin  in  Carderock  Division,  NSWC,  following  standard  test 
procedures  for  barehull  resistance  test,  wake  measurements  in  the  propeller  plane,  propeller 
performance  tests  in  open  water,  and  self-propulsion  tests.  Borda*  reported  the  details  of 
experimental  procedures  and  the  results. 

RESULTS  AND  DISCUSSION 


Figure  4  shows  the  discretized  ship  hull.  320  panels  were  used  to  approximate  one  side  of 
the  ship  hull  surface:  28x8=224  panels  for  main  hull  and  12x8=96  panels  for  bow  dome.  Using  a 
double  body  model  we  examine  how  the  actuator  disk  model  influences  the  local  flow 
characteristics.  Figure  5  shows  the  computed  wake  in  terms  of  cross  flow  velocity  vectors  at  the 
port  side  of  double  body  model  at  Froude  number  0.414.  All  length  scales  are  nondimensionalized 
by  half  ship  length(Lpp/2.0).  The  propeller  radius  Rp  becomes  0.03648  and  the  hub  radius  Rh  is 
0.01028.  TTie  port  side  propeller  is  located  at  P(x,y,z)=  (-0.89178,-0.06,-0.07506).  The  velocities 
are  computed  at  four  different  cross  planes  situated  at  x=  -0.8188,  -0.8553,  -0.8735,  and  -0.8845 
and  these  x-coordinates  are  equivalent  to  2.0* Rp,  1.0*Rp,  0.5  *Rp,  and  0.2*Rn  distance  away 
(upstream)  from  the  actuator  disk.  At  each  cross  plane,  4  different  radii:  r=l.l*Rp,  1.0*Rp, 
0.7*Rp,  and  0.5*Rp  are  considered  and  at  each  radius  the  velocity  components  Vx,Vy,  and  Vz  are 
computed  at  24  points  equally  divided  along  the  circumference.  Figure  5  has  a  pair  of  arrow  plots: 


*  Borda,  Gary  C.,  David  Taylor  Model  Basin,  CDNSWC,  in  a  document  of  higher  classification 
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one  with  and  the  other  without  an  actuator  disk.  As  shown  in  Figure  5,  the  differences  in 
magnitude  and  direction  between  each  pair  become  significant  as  the  local  flow  approaches  the 
actuator  disk. 

In  Table  2,  the  axial  and  total  velocities  are  compared  between  with  and  without  an 
operating  propeller.The  values  of  each  velocity  are  obtained  by  arithmetic  mean  at  each  cross 
section.  Without  operating  propeller  both  axial  and  total  velocities  decrease  their  magnitudes 
progressivley  along  the  x  axis.  An  operating  propeller  accelerates  local  flow  and  the  flow  speed 
increases  as  it  approaches  to  the  propeller.  The  computed  results  show  about  5.5%  increment  in 
speed  at  the  cross  section  x=  -0.8845. 

786  panels  were  used  to  cover  the  free  surface  domain:  58x12=696  panels  for  the  main 
portion  of  the  free  surface  and  15x6=90  panels  for  the  free  surface  following  the  transom  stem.The 
total  1 106  panels  were  used  to  represent  one  side  of  the  plane  of  symmetry. 

Model  experiments  were  conducted  at  Froude  number  ranging  from  0.138  to  0.440  and  the 
wave-making  resistances  were  obtained  from  the  wave  cut  measurement.  In  Figure  6,  the 
computed  wave-making  resistances  are  compared  with  the  measurements.  At  Froude  number 
higher  than  0.276  the  computed  results  show  good  agreement  with  the  measurements,  but  at 
Froude  number  lower  than  0.276  some  difference  is  noticeable.  As  mentioned  in  free  surface  part 
of  numerical  approach,  the  transom  stem  was  wet  until  the  model  reached  the  Froude  number 
0.276.  We  are  interested  in  ship  performance  at  high  speeds,  and  assumed  dry  transom  stem  in 
computation  throughout  the  considered  Froude  number  range. 

In  Table  3,  the  computed  potential  fractions  for  a  double  body  model  are  given  at  four 
Froude  numbers;  0.276,  0.3447,  0.414,  and  0.44.  It  is  well  known  from  D'Alambert's  paradox 
that  the  resistance  of  a  body  in  unbounded  potential  flow  has  to  be  zero.  The  non-zero  resistance 
coefficient  0.10452xl03  without  propeller  comes  from  the  fact  that,  as  discussed  in  the  free 
surface  part  of  numerical  approach,  the  double  body  is  not  closed  completely  at  the  stem  and  left 
open  intentionally  to  accommodate  dry  transom  stem  condition.  The  computed  results  show  that 
the  contribution  of  potential  fraction  to  the  total  thrust  deduction  fraction  is  51%  at  Fn=0.276,  and 
increases  as  speed  increases,  and  becomes  66%  at  F„=0.44. 

Table  4  shows  the  contribution  due  to  the  free  surface  wave  fraction.  About  30% 
contribution  is  observed  at  F„=0.276  and  0.3447, 25%  at  F„=0.414,  and  26%  at  Fn=0.44.  Table  5 
summarizes  the  thrust  deduction  fractions  obtained  from  the  self-propulsion  model  experiments 
and  from  the  computed  potential  and  wave  fractions.  The  experimental  t  values  include  all  three 
components;  potential(tp),  viscous(tv),  and  free-surface  wave(tw),  but  the  computations  are  limited 
to  estimates  of  tp  and  tw.  The  following  observations  can  be  made: 

(1) At  the  Fn's  between  0.276  and  0.440,  the  computed  tp+tw  values  lie  between  80%  and 
92%  of  the  measured  thrust  deduction.  The  difference  may  be  caused  by  the  neglected 
viscosity  effect,  or  by  the  limits  of  computational  accuracy  and  capability,  or  both. 

(2)  The  component  tw  has  been  customarily  neglected  because  of  the  numerical  difficulties 
involved  in  handling  the  free  surface  condition  and/or  because  it  was  thought  small 
quantity.  The  computational  results  indicate  that  for  this  particular  model  the  tw  term  is  not 
small,  between  25  and  30%  of  the  measured  t,  and  should  not  neglect. 

(3)  The  measured  wakes  lie  between  0.01  and  0.04  at  the  considered  speeds.  These  values 
are  extremely  small  compared  with  those  values  of  0.2  and  0.5  which  are  common  in 
typical  low  and  high  block  coefficient  commercial  ships,  respectively.  For  this  ship  fitted 
with  the  inclined  propeller  shaft  with  respect  to  the  keel  line  the  measured  values  indicate 
that  the  propeller  plane  seems  completely  outside  of  the  developed  hull  boundary  layer. 
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(4)  The  thrust  deduction  fractions  are  also  sensitive  to  hull  forms.  The  measured  t  values  for  a 
cruiser  stem  type  ship  reported  in  [2]  and  [5]  varies  between  0. 1  and  0.3  while  those  are 
between  0.07  and  0.095  for  this  transom  stem  ship. 

CONCLUSIONS  AND  RECOMMENDATIONS 

From  the  comparisons  between  the  experimental  and  the  computed  results,  it  has  been 
demonstrated  that  the  panel  method  is  capable  of  handling  an  important  practical  hull  configuration 
and  can  provide  reasonable  qualitative  and  quantitative  results  for  the  potential  and  free-surface 
wave  components  of  thrust  deduction.  Thus  far  the  present  panel  method  approach  for  the 
calculation  of  the  free  surface  effect  on  thrust  deduction  has  been  applied  in  detail  to  this  one 
example  case.  Though  the  viscous  effects  are  not  considered  here,  a  method  is  now  in  hand  to 
make  estimate  of  the  effect  of  the  free  surface  waves  on  the  thrust  deduction  fraction.  With  regard 
to  these  quantitative  results,  the  following  general  comments  apply: 

-  Over  the  speed  range  of  the  ship,  the  free  surface  wave  component  of  thrust  deduction,  tw , 
is  not  negligibly  small. 

-  For  a  transom  stem  ship  with  some  inclined  propeller  shaft  angle,  the  measured  onset  flow 
to  the  propeller  seems  essentially  uniform.  The  propeller  disk  is  lying  outside  the  hull 
wake  shadow  and  the  viscous  fraction  tv  is  expected  small. 

-  To  improve  the  prediction  of  the  local  propeller-induced  flow  field,  the  simple  actuator  disk 
model  with  constant  sink  strength  must  be  replaced  with  a  more  accurate  propeller  model 
which  could  include  the  effect  of  onset  flow  inclination  as  well. 


ACKNOWLEDGMENTS 

The  author  expresses  his  gratitude  to  Dr.  M.B.  Wilson  and  Dr.  K.H.  Kim  for  their  valu¬ 
able  suggestions  and  support 


10 


Figure  2.  Typical  Panel  Arrangement  for  a  Ship  and  Free  Surface 


Figure  4.  Ship  Hull  Panel  Network 


(w/o  prop,  x=0.5R) 

PLOT  OF  THE  SECONDARY  FLOW 

(with  prop,  x=0.5R) 
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)  tiS  /  f 

run  auu  f A,  o* 
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PLOT  OF  THE  SECONDARY  FLOW 

PLOT  OF  THE  SECONDARY  FLOW 
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run  khz  ?  a,  rr 

Figure  5  -  continued 


Figure  6.  Comparison  of  Wave-Making  Resistance  Coefficients  Predicted  by  SWIFT  and 
from  Model  Experiments 
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Table  1  Ship  Particulars 


B(beam)/L(hull  length) 


D(draft)/B 


Dome  Length/L 


Block  Coeff. 


Wetted  Surface  Area  Coeff. 


Waterplane  Area  Coeff. 


Transom  Stem  Area  Coeff. 


where  Block  Coeff.  =  Volume/(L*B*D) 

Wetted  Surface  Area  Coeff.  =  Area/(L*(B+2D» 
Waterpalne  Area  Coeff.  =  Area/(L*B) 

Transom  Stem  Area  Coeff.  =  Area/(B*D) 


Table  2.  Mean  Axial  and  Total  Velocity  Increment  due  to  Actuator  Disk 

at  F„  «  0.0.4140 


Section  I 

Without 

’ropeller 

j  With  Propeller  i 

Velocity 

ncrement(%! 

i 

tt  X. 

|Vx| 

M 

|Vx| 

M 

A|Vx| 

A|V| 

-0 

.8188 

0.97909 

0.98200 

0.99213 

0.99542 

1.33 

1.37 

-0.8553 

0.97459 

0.97694 

0.99801 

1.00088 

2.40 

2.45 

-o 

.8735 

0.97257 

0.97464 

1.00954 

1.01227 

3.80 

3.86 

-0.8845 

0.97144 

0.97332 

1.02471 

1.02778 

5.48 

5.60 

where  IVh 

=  sqrt(  Vx2  + 

Vy2  +  Vz2)  and  the  propeller  is  located  at  x  = -0.89178. 
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Tabic  5.  Thrunst  Deduction  Comparison  between  Experiment  and  Computation 

Thrust  Deduction  Fraction 

Fn 

Gjh 

*e 

Computation 

Ratio  (%) 

(Exp) 

(Exp) 

*P 

t* 

tp  +t* 

t 

(tp+tw)/t 

0.2760 

0.3186 

0.03 

0.04877 

0.02837 

mm 

0.095 

81.20 

0.3447 

0.3466 

0.04 

0.04895 

0.02885 

mm 

0.095 

81.90 

0.4140 

0.4828 

0.025 

0.04713 

0.02099 

0.06812 

0.085 

80.14 

0.4400 

0.5395 

0.01 

0.04610 

0.01829 

HB9 

0.070 

91.99 
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